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We study the influence of Gribov copies on gluon and ghost propagators in lattice Landau gauge. For the gluon 
propagator, Gribov noise seems to be of the order of magnitude of the numerical accuracy. On the contrary, for 
the ghost propagator, Gribov noise is clearly observable, at least in the strong-coupling regime. We also observe, 
in the limit of large lattice volume, a gluon propagator decreasing as the momentum goes to zero. Finally, we 
introduce an implementation of the method of Fourier Acceleration which avoids the use of the fast Fourier 
transform, being well suited for parallel and vector machines. We apply it to the case of Landau gauge fixing, 
and study its performance on APE computers. 



1. GRIBOV NOISE 

In ref. (l]] Gribov showed that, for non-abelian 
gauge theory, the standard gauge-fixing condi- 
tions used for perturbative calculations do not 
fix the gauge fields uniquely. The existence of 
these Gribov copies does not affect the results 
from perturbation theory, but their elimination 
could play a crucial role for non-perturbative fea- 
tures of these theories. 

In lattice gauge theories gauge fixing is, in prin- 
ciple, not required. However, because of asymp- 
totic freedom, the continuum limit is the weak- 
coupling limit, and a weak-coupling expansion re- 
quires gauge fixing. Thus, one is led to consider 
gauge-dependent quantities on the lattice as well. 
Unfortunately gauge fixing on the lattice is af- 
flicted by the same problem of Gribov copies en- 
countered in the continuum case Q . 

In order to get rid of Gribov copies the physical 
configuration space has to be identified with the 
so-called fundamental modular region A, which is 
defined (in the continuum) as the set of absolute 
minima of the functional [|3| 
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Similarly, on the lattice, we can eliminate Gribov 
copies looking for the absolute minimum of the 
functional £jj[g] (minimal Landau gauge) B 
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'Poster presented by A. Cucchieri. 



Given the appearance of Gribov copies in nu- 
merical studies, we need to understand their in- 
fluence (Gribov noise) on the evaluation of gauge- 
dependent quantities. To this end we compare 
the results for gluon and ghost propagators us- 
ing two different averages §: the average con- 
sidering only the absolute minima (denoted by 
"am"), which should give us the result in the min- 
imal Landau gauge; and the average considering 
only the first gauge-fixed gauge copy generated 
for each configuration (denoted by "fc"). The lat- 
ter average is the result that we would obtain if 
Gribov noise were not considered. 

Our data [|[ Table 2] show absence of Gribov 
noise for the gluon propagator. In fact, data cor- 
responding to the minimal Landau gauge (abso- 
lute minima) are in complete agreement, within 
statistical errors, with those obtained in a generic 
Landau gauge (average "fc"). This happens even 
at /? = 0, where the number of Gribov copies is 
very large and Gribov noise, if present, is more 
easily detectable. On the contrary, a nonzero 
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Figure 1. Plot of the three- momentum-space gluon 
propagator D(k) ("fc"-data), as a function of the 
square of the lattice momentum p 2 (k). Data corre- 
spond to V = 16 3 (□), V = 32 3 (O) and V = 64 3 
(*), at (3 = 5.0. Error bars are one standard devia- 
tion. Averages are taken over 40 gauge-fixed config- 
urations. 



Gribov noise for the ghost propagator can be 
clearly observed [|[ Table 3]. In particular, data 
corresponding to the absolute minima (average 
"am") are constantly smaller than or equal to 
the corresponding "fc"-data. This effect is small 
but clearly detectable for the values of f3 in the 
strong-coupling region. (This was not observed 
at (3 = 2.7. However, at this value of (3 almost no 
Gribov copies were produced, even for a lattice 
volume V = 16 4 , and therefore we cannot expect 
a difference between the two sets of data.) This 
result can been qualitatively explained ^| . As for 
the infrared behavior of these two propagators, 
the data for the ghost propagator show a pole 
"between" the zeroth-order perturbative behav- 
ior p~ 2 — valid at large momenta — and the 
singularity predicted in , but in agreement with 
the pole p _2 ' 1+s ^(0 < s < 1) recently obtained in 
ref. ||. For the gluon propagator the data show, 
in the strong-coupling regime, a propagator de- 
creasing as the momentum goes to zero. This 
anomalous behavior, predicted in jjJHQl, is still 
observable at (3 = 1.6, if large volumes are consid- 
ered [p| . This result is also observable in the scal- 



ing region in the three-dimensional case, and in 
the limit of large lattice volume (see Fig. Q). Fi- 
nally, the behavior of the zero three-momentum- 
space gluon propagator is strongly affected by the 
zero-momentum modes of the gluon held ^ Fig. 
2], as predicted in jjj). 

2. FOURIER ACCELERATION 

In order to minimize the functional £jj[g\ de- 
fined in eq. (||), and to reduce critical slowing- 
down, we can use the Fourier accelerated algo- 
rithm Jl^,|ll| • With this algorithm the update is 
given by g^ ew \x) = R{x) g iold) (x), where 
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Here F is a Fourier transform, a is a tuning pa- 
rameter, p 2 (k) is the square of the lattice mo- 
mentum, and V • A is the lattice divergence of the 
gluon held A^. However, this algorithm is of dif- 
ficult implementation on parallel machines, due 
to the use of the fast Fourier transform (FFT). 

Let us notice that F~ 1 p- 2 (k)F = (-A) -1 , 
where A is the lattice Laplacian operator. Thus, 
the FFT can be avoided by inverting A using an 
algorithm that requires the same computational 
work (i.e. VlogN), such as a multigrid (MG) al- 
gorithm with W cycle and piecewise-constant in- 
terpolation. At the same time, using MG, we can 
reduce the computational work with a good initial 
guess for the solution, and we can choose the ac- 
curacy of the solution. (With FFT the accuracy 
is hxed by the precision used in the numerical 
code.) We note that the tuning parameter a is 
usually hxed with an accuracy of a few percent, 
and thus the inversion of A should not require a 
very high accuracy either. 

We started our simulations on an IBM RS- 
6000/340 workstation. We tested different types 
of multigrid cycles: 7 = (Gauss-Seidel update), 
7 = 1 ( V cycle) and 7 = 2 (W cycle) . We see that 
MG with 7 = 2, two relaxation sweeps on each 
grid, a minimum of two full multigrid sweeps for 
each version of A, and an accuracy of 10~ 6 , is 
equivalent to an FFT algorithm jig]. 

In Table | we report results obtained at (3 = 00 
for the FFT algorithm, and for MG. Clearly, the 
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algorithm 


V 


GF-sweeps 


CPU-time 


FFT 


4 4 


13.6 ±0.3 


36 ±1 


MG 


4 4 


13.7 ±0.3 


65 ±3 


FFT 


8 4 


16.1 ±0.3 


972 ± 16 


MG 


8 4 


16.1 ±0.3 


1381 ±49 


FFT 


16 4 


20.0 ±0.4 


25254 ±560 


MG 


16 4 


19.9 ±0.4 


33797 ± 1114 



Table 1 : Comparison between the FFT and MG al- 
gorithms at (3 = oo. 



two algorithms have a similar performance, show- 
ing a number of gauge-hxing sweeps increasing 
logarithmically with the lattice size N, and the 
CPU-time increasing as N 4 \ogN. We also did 
a test at (3 = 2.2 (see Table |). Again, MG is 
equivalent to the FFT algorithm. 



algorithm 


V 


GF-sweeps 


CPU-time 


FFT 


8 4 


333 ± 27 


19998 ± 1236 


MG 


8 4 


312 ±25 


18722 ± 1496 



Table 2 : Comparison between the FFT and MG al- 
gorithms at (3 = 2.2. 



In order to parallelize, the idea is to use as the 
coarsest grid for the multigrid algorithm a grid 
with volume equal to or larger than the number 
of nodes of the parallel machine. For example, 
for an APE100 computer with 8 3 nodes we im- 
plemented MG with the coarsest grid 8 4 . Then, 
on the coarsest grid, we can use a Gauss-Seidel re- 
laxation if its volume is small. Otherwise we can 
use a Conjugate Gradient algorithm to relax the 
solution. (This combination MG±CG has been 
used in the past to accelerate MG on vector ma- 
chines p3|.) In this way, the computational work 
for the inversion of A still increases as VlogN, pro- 
vided that we keep fixed the size of the coarsest 
grid. We tested this combination first on a work- 
station for a 16 4 lattice with coarsest grid 8 4 at 
(3 = oo, performing two CG-sweeps when relaxing 
on the coarsest grid. We obtained 20.0 ± 0.4 for 
the GF-sweeps, and 46281 ± 1732 for the CPU- 
time. Therefore (see Table |l|) the performance of 
this gauge-fixing algorithm is essentially equiva- 
lent to that of FFT and MG. 

Similarly, the performance at (3 — 2.2, for an 8 4 
lattice with coarsest grid 4 4 , is comparable to the 



performance of the FFT and MG algorithms (see 
Table ||): we obtained 314±25 for the GF-sweeps, 
and 23599 ± 1884 for the CPU-time. 

Finally, we implemented the MG±CG algo- 
rithm on an APE 100 computer comparing its per- 
formance with a standard overrelaxation (OVE) 
and an unaccelerated local algorithm (the so- 
called Los Alamos algorithm, LOS) (TTJ . The 
number of gauge-fixing sweeps obtained, at /3 = 
oo and for lattice volume 16 4 , was 131±3 for LOS, 
34.8 ± 0.5 for OVE, and 16.4 ± 0.1 for MG+CG. 
Clearly the MG±CG algorithm is able to reduce 
the number of gauge-fixing sweeps compared to 
the two local algorithms. 

We plan to extend the tests on APE computers 
to larger lattice volumes. 
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